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ABSTRACT 

The sdB star KPD 0422+5421 was discovered to be a single-lined spectroscopic binary 
with a period of P = 0.0901795 ±(3 x 10" 7 ) days (2 hours, 10 minutes). The U and B 
light curves display an ellipsoidal modulation with amplitudes of sa 0.02 magnitudes. 
The sdB star contributes nearly all of the observed flux. This and the absence of 
any reflection effect suggest that the unseen companion star is small (i.e. i? C omp ~ 
0.01 Rq) and therefore degenerate. We modeled the U and B light curves and derived 
i = 78.05 ±0.50° and a mass ratio of q = M comp /M sdB = 0.87 ±0.15. The sdB star fills 
69% of its Roche lobe. These quantities may be combined with the mass function of the 
companion (f(M) = 0. 126 ± 0.028 M ) to derive M sdB = 0.72 ± 0.26 M© and M comp = 
0.62 ± 0.18 Mq. We used model spectra to derive the effective temperature, surface 
gravity, and helium abundance of the sdB star. We found T c g — 25, 000 ± 1500 K, 
log<? = 5.4 ± 0.1, and [He/H] = —1.0. With a period of 2 hours and 10 minutes, 
KPD 0422+5421 has one of the shortest known orbital periods of a detached binary. 
This system is also one of only a few known binaries which contain a subdwarf B 
star and a white dwarf. Thus KPD 0422+5421 represents a relatively unobserved, and 
short-lived, stage of binary star evolution. 

Key words: binaries: close — stars: variable — stars: individual (KPD 0422+5421) 



1 INTRODUCTION 

The star KPD 0422+5421 was discovered by Downes (1986) 
during the course of a search for very blue stars in the galac- 
tic plane, and classified as an sdB star, i.e. a hot, hydrogen- 
rich subdwarf. R. Saffer (private communication) derived 
from spectroscopic measurements an effective temperature 
T e g = 26,050 and a surface gravity \ogg — 5.51; Saffer did 
not communicate a helium abundance. The sdB classifica- 
tion prompted one of the current authors to include KPD 
0422+5421 in a list of candidates to be monitored for rapid 
pulsations (see e.g. Kilkenny 1997a). During the course of 
the mandatory 90 minute high speed photometric run on 
the star, it appeared to vary with a period of about an hour, 
at an amplitude of roughly 0.01 mag. Further spectroscopic 
and photometric observations show that KPD 0422+5421 
is a close binary star with a period of 2.1643 hours, one of 
the shortest known orbital periods among the detached bi- 
naries. We argue that the companion star is a white dwarf, 
making KPD 0422+5421 one of a small number of known 
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systems with a sdB and a white dwarf. We report here our 
observations, data reductions, and data analysis. 



2 THE HIGH SPEED PHOTOMETRIC 
OBSERVATIONS AND ANALYSIS 

All the photometry reported here was obtained at the 
University of Texas' McDonald Observatory on Mt Locke, 
Texas. Two different instrumental configurations were 
used: the Louisiana State University Photometer (P-LSU) 
mounted on the 0.9m telescope, and the Stiening Photome- 
ter attached to the 2.1m telescope. The P-LSU is a stan- 
dard two-channel photoelectric photometer (see e.g. Grauer 
& Bond 1981), equipped with a blue-sensitive Hamamatsu 
R647 photomultiplier tube. The Stiening Photometer is 
a four-channel instrument which allows simultaneous high 
speed photometry in four different wavebands (Home & 
Stiening 1985; Wood, Zhang & Robinson 1993). The cur- 
rent passbands of the instrument are similar to Johnson's 
UBVR, although there are notable differences (Wood et al. 
1993). One channel of a second photometer (P45) was used 
as a fifth channel with the Stiening photometer, to monitor 



© 1998 RAS 



2 Koen et al. 



Table 1. Log of the observations, all of which were made at the 
University of Texas' McDonald Observatory. All 0.9m observa- 
tions were obtained under good photometric conditions, while all 
observations on the 2.1m telescope were acquired under poor con- 
ditions. The instrumental configurations are described in the text. 
JD 2,450,752.8437 is 1997 October 31.3438 UT 



Starting Time Run Length Telescope 



JD 2,450,000+ 


(Hours) 


(meters) 


752.8430 


1.4 


0.9 


753.8599 


3.0 


0.9 


754.7085 


6.9 


0.9 


785.7645 


4.0 


2.1 


787.7779 


3.9 


2.1 



Table 2. A listing of the largest peaks in the amplitude spectrum 
of all the high speed photometric data. 



Frequency 


Period 


Amplitude 


(c/d) 


(days) 


(mmag) 


22.055 


0.045341 


14.4 


22.086 


0.045277 


15.7 


22.117 


0.045215 


16.6 


22.148 


0.045151 


17.0 


22.178 


0.045088 


17.0 


22.209 


0.045027 


16.6 


22.240 


0.044964 


15.8 


22.271 


0.044901 


14.5 
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Time (hours since HJD 3,450.754.708) 
Figure 1. The white light observations obtained on JD 2,450,754 
(1997 November 2 UT) using the 0.9m telescope. One minute 
averages of the ten-second integrations are shown. The zeropoint 
has been set to the data mean value. 



a comparison star. Observations with both P-LSU and P45 
photometers were in "white light", i.e. no niters were placed 
in the light beam. The resulting effective wavelengths are 
similar to Johnson's B, but with substantially wider band- 
passes. Integration of 10 and 5 seconds were used when run- 
ning respectively the P-LSU and Stiening photometers. 

The photometric runs are cataloged in Table 1; approxi- 
mately 11 and 8 hours of data were obtained on the 0.9m and 
2.1m telescopes respectively. The longest continuous light 
curve is plotted in Figure lL There is a clear modulation 
with a period of « 65 minutes. Amplitude spectra of vari- 
ous combinations of photometric runs are shown in Figure 
|i| plotted over the frequency range of greatest interest. The 
low frequency content of the amplitude spectra has been 
removed to some extent by detrending the observations lin- 
early. In addition, the fifth channel observations have been 
used to correct the data acquired on JD 2450785 for large 
atmospheric transparency drifts. 

The one cycle per day aliasing pattern is clearly vis- 
ible in the top panel of Figure ^, which shows the am- 
plitude spectrum of the combined 0.9m data (from three 
successive nights). The aliasing pattern is much worse for 
the combined two 2.1m runs (middle panel), since these 
were short. The bottom panel shows the amplitude spec- 
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Figure 2. The amplitude spectrum of the all the 0.9m data (top 
panel); all the 2.1m data (middle panel); and all the data from 
both telescopes (bottom panel). The observations from each night 
were individually detrended by a first order polynomial. In the 
case of the 2.1m observations, fluxes from the four channels have 
been co-added. 



trum of all the data; although not visible on the scale of 
Figure ^, the main peak is composed of an alias pattern 
of sub-peaks, spaced 0.031 cycles d" 1 apart. The frequen- 
cies for the peaks with the largest amplitudes are listed in 
Table 2. A nonlinear least squares fitting procedure, with 
a starting guess for the frequency corresponding to the 
largest-amplitude alias gave / max = 22.178 ± 0.00016 cycles 
day" 1 [P = 0.0450897 ± (3 x 10-7) days] and amplitude 
= 17.2 ± 0.3 mmag. In spite of the small formal frequency 
error, we would hesitate to choose the correct alias based on 
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Table 3. Amplitudes and phases (with respect to the first obser- 
vation) of the variations in the different wavebands. The contents 
of this Table were derived by fitting a sinusoid with / = 22.178 
c/d to the U, B, V and R, data by linear least squares. The formal 
standard errors are given in brackets below each estimate. 



Table 4. The orbital parameters for KPD 0422+5421. 





V 


B 


V 


R 


Amplitude (mmag) 


19.4 


17.3 


14.8 


11.8 


(s.e.) 


(0.3) 


(0.2) 


(0.3) 


(0.3) 


Phase 


1.64 


1.61 


1.65 


1.51 


(s.e.) 


(0.02) 


(0.01) 


(0.02) 


(0.03) 



the photometry alone. Some of the reasons for our misgiv- 
ings are: the two large peaks in the amplitude spectrum have 
almost identical heights; corrections for drifts in the pho- 
tometry zeropoints are uncertain; and the 2.1m data used 
in the procedures described above was the simple sum of the 
fluxes in all four filters. This choice was based on the desire 
to maximise the available flux, but implies a distinct mis- 
match of the wavebands of the two sets of observations. It is 
quite conceivable that the relative heights of the amplitude 
spectrum aliases might change if a different combination of 
Stiening photometer fluxes is used. However, we have ad- 
ditional constraints imposed by spectroscopic observations. 
We show below that the orbital period of KPD 0422+5421 
is P = 0.0901795 ± (3 x 10~ 7 ) days, which corresponds to a 
period of P = 2// max . 

A linear least squares programme was used to fit the 
frequency as determined above to the combined two 2.1m 
runs; this was done separately for each of the four Stiening 
photometer wavebands. The amplitudes are given in Table 3. 
Although formal standard errors of the fits are given, these 
should be viewed only as very rough guides. In particular, 
errors in different waveband datasets are probably highly 
correlated, so that the evaluation of e.g. colour phase differ- 
ences on the basis of the formal errors may be misleading. 



3 SPECTROSCOPIC OBSERVATIONS AND 
ANALYSIS 

We obtained 14 spectra of KPD 0422+5421 starting at HJD 
2,450,756.8586 (1997 November 4.3583 UT) with the 2.7m 
telescope at the McDonald Observatory using the Large 
Cass Spectrograph, a 600^/mm grating (blazed at 4200A), 
and the Til 800 x 800 CCD. The resulting spectral resolu- 
tion is ~ 3.5 A (FWHM) with wavelength coverage of 3525- 
4940 A. The exposure times ranged from 10 to 20 minutes, 
and the signal-to-noise ratios in the final reduced spectra 
were ~ 20 — 40 in the continuum near H/3. There were cir- 
rus clouds present, so the spectra could not be placed on 
an absolute flux scale. Instead, we used nine observations 
of five different spectrophotometric flux standards from the 
nights of November 1-3, 1997 (which were photometric) to 
calibrate out the instrumental response. Since the slit could 
not be rotated to follow the parallactic angle, most of our 
program objects were observed at hour angles near zero. 
The relative flux scale is reasonably accurate between H/3 
and the Balmer jump — however the relative calibration at 



Parameter 



Value 



spectroscopic period (days) 
2 X photometric period (days) 
T (spect) (HJD 2,450,000+) 
T (photo) (HJD 2,450,000+) 
-KsdB velocity (km s ) 
7 velocity (km s -1 ) 
f(M) (M ) 



0.0907 ± 0.0020 
0.0901795 ± (3 x 10~ 7 ) 
756.9381 ± 0.0010 
785.8199 ± 0.0050 
237 ± 18 
-57 + 12 
0.126 + 0.029 



either end of the covered spectral range is suspect in part 
because some of the KPD 0422+5421 spectra were obtained 
at relatively large airmass (1.5-2.0). 

We used the cross-correlation technique of Tonry & 
Davis (1979) to determine the radial velocities of the KPD 
0422+5421 spectra. The cross-correlation functions (CCFs) 
were computed over the wavelength interval 3800-4920 A. 
We used a synthetic spectrum generated from a model atmo- 
sphere (see below) as the template spectrum (although our 
results were almost identical when the first observation was 
used as the template). In each case the CCFs were strong 
and well-defined. The velocities corresponding to the cen- 
troid of the CCF peaks were determined by a parabolic fit 
to the six pixels surrounding the maximum. 

Large radial velocity variations were evident. We fitted 
a four-parameter sinusoid to the 14 radial velocities. The 
best-fitting sinusoid had xt = 6.5, an indication that the 
errors on the radial velocities were probably too small. We 
scaled the errors on the individual velocities to give xt — 1 
for the fit. The resulting spectroscopic elements are given 
in Table 4. Note that the spectroscopic period of P sp cct = 
0.0907+0.0020 days is roughly twice the photometric period 
derived above. The spectroscopic period and the velocity 
semi-amplitude of A" s( jb = 237 ± 18 km s _1 imply a mass 
function of the (unseen) companion star of 



f(M) 



PKt 



M c 3 omp sin 3 i 



2-KG (Afcomp + Af sdB ) 



0.126+0.029 M Q (1) 



where M a 



is the mass of the unseen companion star, 



M s dB is the mass of the sdB star, and % is the orbital inclina- 
tion. The above equation implies M com p > 0.126+0.029 Mq. 

The fact that the photometric period is half the spec- 
troscopic period suggests that the optical modulations we 
observe are due to the well-known "ellipsoidal" variations. 
The cause of the ellipsoidal modulations is easy to under- 
stand. The sdB star we observe is slightly distorted by the 
unseen companion star. As the sdB star moves around in its 
orbit, its projected area on the sky changes. The changes in 
the projected area give rise to the observed changes in the 
total flux. We expect to observe the maximum flux at the 
quadrature phases when the sdB star is viewed "side-on." 
Similarly, we will see minima in the light curve during the 
conjunction phases when the sdB star is viewed "end-on." 
The photometric minimum at the spectroscopic phase 0.25 
(the superior conjunction of the sdB star) will be slightly 
deeper due to gravity darkening (see Avni 1978 and the 
discussion below). Thus we can use the spectroscopically 
determined phase to select the correct alias period from 
among the several of the most likely candidates. We find 
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Figure 3. Top: The folded radial velocities of KPD 0422+5421 
and the best-fitting sinusoid. Middle: The folded light curve from 
the 0.9m (white light), binned into 50 phase bins. Bottom: The 
folded light curve from the 2.1m (B band), binned into 100 phase 
bins. The errors shown for the binned light curves are the errors 
of the mean in each bin. 



that only the period of 2// max = 0.0901795 ± (3 x 10~ 7 ) 
days correctly phases the 0.9 and 2.1 meter photometric 
data with respect to the radial velocity curve. The time of 
the superior conjunction of the sdB star, as determined from 
the 2.1m data, is T (photo) = HJD 2, 450, 785.8199 ± 0.005. 
The values of To (photo) and Tb(spect) differ in phase by 

= 320.270 ± 0.013, consistent with the expected value 
of A(/> = 320.25. The folded radial velocity curve and the 
folded light curves from the 0.9m (white light) and the 2.1m 
(B) are shown in Figure |§|. Models of the ellipsoidal light 
curve will be discussed in Section ^ below. 

We created a "restframe" spectrum of KPD 0422+5421 
by Doppler correcting all 14 spectra to zero velocity and av- 
eraging them. This restframe spectrum is displayed in Figure 
^| Two important points should be noted. First, the spec- 
trum does not show any obvious lines from the companion, 
quite unlike typical "composite" spectrum subdwarf stars 
from the Palomar-Green (PG) survey (Ferguson, Green, & 
Liebert 1984; Green, Schmidt, & Liebert 1986; Orosz, Wade, 
& Harlow 1997). Second, the He lines are relatively weak. 
We will discuss models of the spectrum in Section below. 



4 CONSTRAINTS ON THE SYSTEM 
GEOMETRY FROM LIGHT CURVE 
MODELS 

We analyzed the U and B light curves from the observations 
on the 2.1m telescope, by using the Wilson- Devinney (W-D) 
code (Wilson & Devinney 1971; Wilson 1979). Only the U 
and B light curves were used in the modelling process as 
these were less noisy than the V and R data. This may be 
due in part to the greater difficulty in correcting the obser- 
vations through the two red filters for random atmospheric 
transparancy changes. 

We can draw some basic conclusions about the system 
even before we begin the modelling. The mass of the sdB 
star is probably not very different from 0.5 Mq (Saffer et al. 
1994). We can compute the orbital separation as a function 
of the total system mass from Kepler's third law: 



3 P 2 GM tota i 
a = • 

We find for KPD 0422+5421 that 
» 0.846 ('%^ X 



(2) 



(3) 



For mass ratios near 1, the radius of the companion's Roche 
lobe is approximately 38% of the orbital separation (Eggle- 
ton 1983), which for a total system mass near 1 Mq, is about 
0.35 R Q for KPD 0422+5421. If the companion star were on 
the main sequence, then it would be later type than about 
M2 in order to have a radius small enough to fit inside the 
Roche lobe (Gray 1992). Even then, it would still fill a siz- 
able fraction of its Roche lobe. Hence we would expect to 
see a reflection effect where the irradiation of the M star by 
the hot sdB star causes extra light to be observed near the 
spectroscopic phase 0.75, as in the case of HW Vir (Wood 
et al. 1993) . It is clear from Figure |^ that the light curves 
are ellipsoidal — there is no hint of a reflection effect. This 
rules out the possibility that the companion is a cool main- 
sequence star. The mass function rules out a brown dwarf 
of smaller radius. In principle, there would not be any siz- 
able reflection effect if the two stars had nearly the same 
temperature. However, if the companion star had a radius 
and temperature similar to the sdB star, we would observe a 



double-lined spectroscopic binary since L c 



The 



system is single-lined, however. We conclude that the com- 
panion star is much smaller than the sdB star, and therefore 
probably is a white dwarf. 

The light curves were fit using "mode 2" of the W-D 
code where "star 1" was the sdB star. Mode 2 is the normal 
mode used when neither star is in contact with its Roche 
lobe. We assumed circular orbits and synchronous rotation, 
both reasonable assumptions considering the proximity of 
the two stars. The monochromatic fluxes were computed as- 
suming blackbodies. We fixed the polar temperature of the 
sdB star at T\ = 26000 K, and assigned limb darkening 
coefficients of u\ = 0.31 and 0.29 for the U and B bands, 
respectively (Wade & Rucinski 1985). The sdB star has a 
radiative envelope, so its gravity darkening exponent g\ was 
set to 1 and its bolometric albedo A\ for reflective heating 
and re-radiation was set to 1. Since the white dwarf con- 
tributes less than 0.1% of the total flux, the values of T%, 
U2, 52, and A2 are relatively unimportant. Thus we allowed 
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Figure 4. The rest frame spectrum of KPD 0422+5421, which is 
the average of the 14 spectra Doppler corrected to zero velocity. 
The signal-to-noise ratio is about 130 near H/3. 



T2 to be a free parameter and set 112 = lii, g% = gi, and 
A2 = Ai. The free parameters in the model fits are the in- 
clination i, the mass ratio q = M CO m P /M s dB, T2, and the 
surface potentials Qi and f?2. The effective wavelengths of 
the Stiening U and B niters were taken to be 3460 A and 
4370 A, respectively (Wood et al. 1993). 

We used a variation of the "gridls" program given in 
Bevington (1969) to optimize the parameter values. We de- 
fined the optimal light curve solution to be the one which 
minimized the \ 2 fit to both the normalized U and B curves 
(in magnitude units) simultaneously. We also used the "dif- 
ferential corrections" routine of the W-D code to check our 
results. We are confident we have found the global \ 2 min- 
imum (rather than a local minimum) since we started the 
optimization routine from a large number of widely sepa- 
rated regions in parameter space. We show in Figure ^ the 
best fit solution to the U and B light curves and the resid- 
uals. The overall fit is reasonably good with Xmin = 364.36 
for 200 data points. The residuals in the U and B filters 
have standard deviations of 0.0029 and 0.0025 magnitudes, 
respectively. Table 5 summarizes the assumed input param- 
eters and the derived parameters from the model light curve 
fits. 

The statistical errors of the parameters derived from the 
light curve models were derived in the following way. The 
parameter of interest (call it Oi) was slightly altered from its 
optimal value (a} — a® + S) and held fixed while the other 
free parameters were adjusted to give the minimum x 2 - The 
process was repeated (a 2 = a\ + 5, etc.) until \ 2 differed 
substantially from Xmin- Figure ^ shows the \ 2 vs - * an d the 
X 2 vs. q curves. A change in the inclination z of ~ 0.5° is 
required to force x 2 to change by 1. We adopt cr, = 0.50°, 
although we note that this is a rough approximation since 
the curve is not parabolic and Xmin i s much larger than the 
number of data points fit. Based on the x 2 vs - 1 curve, we 
adopt a q = 0.15 for the sake of discussion (and with the same 
caveats discussed for <7j). The quoted errors on the other 
parameters given in Table 5 were derived assuming Gaussian 
errors. Eggleton's (1983) formula was used to compute the 




0.5 1 1.5 2° 



Photometric Phase 



Figure 5. From top to bottom: The U light curve and the best 
fit model; the U residuals in the sense of data minus model; the 
B light curve and fit; and the B residuals. Each point has been 
plotted twice. See Table 5 for the light curve model parameters. 



effective radii of the Roche lobes 

RmjsdB) 0.49g- 2/3 

a ~ 0.6g~ 2 / 3 +ln(l + q- 1 /3) 

fla(comp) = 0.49g 2/3 

a 0.6g 2 /3 + ln(l + q 1 / 3 ) ' [ ' 

The W-D code computes the radii of the components using 
the values of Qi, Q2, and the orbital separation a. The sdB 
star is well within its Roche lobe. 

We tried model fits with the temperature of the sdB star 
set to 25,000 K and to 27,000 K. The best-fit values of i and q 
were nearly identical to those values given in Table 5. In that 
same vein, we also tried model fits where the limb darkening 
coefficients were set to 0.25 and 0.35 for both filters. Again, 
the best-fit values of i and q did not change significantly. 
Finally, we found that the best-fit values of i and q did not 
change significantly when the effective wavelengths of the 
U and B filters were changed by 100 A. We conclude that 
our results presented in Table 5 are insensitive to the exact 
values chosen for the sdB parameters. 

We emphasize that we have only estimated the statisti- 
cal errors and not any systematic errors. There are no doubt 
systematic errors present, given the small light curve ampli- 
tudes and the nature of the observations (i.e. high speed 
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Figure 6. Left: The x 2 vs - i curve. The value of i was fixed at each 
value along the curve and the other parameters were adjusted 
as to give the lowest \ 2 - The dashed-dotted lines denote Xmim 
Xmin + 1, Xmin + 4, and Xmin + 9. Right: The x 2 vs. q curve. 



Table 5. Fitted parameters for KPD 0422+5421. 



Parameter Value 



i (deg) 78.02 ± 0.50 

g 0.87 ±0.15 

T 2 (K) ss 4040 

ill 4.8642 ±0.24 

tt 2 « 57.6086 

M sdB (M Q ) 0.72 ±0.26 

Mcomp (Mq) 0.62 ±0.18 

-RsdB {Re) 0.248 ±0.008 

R corap (Rq) « 0.01 

logg sd B (cgs) 5.50 ±0.21 

loggcomp (cgs) fa 7.9 

R m (sdB) (Rq) 0.36 ±0.05 

K fll (comp) (Rq) 0.35 ±0.03 

expected V TO t s'mi (km s — 1 ) 137 ± 4 

a {Rq) 0.93 ±0.11 



photometry). However, the white light curve from the 0.9m 
and the B light curve from the 2.1m are very similar in am- 
plitude, and their relative phases are as expected. There is 
also nothing unusual about the derived component masses. 
The mass of the sdB star (0.72 ±0.26 Mq) is consistent with 
the "canonical" extended horizontal branch mass of 0.5 Mq 
(Caloi 1989; Dorman, Rood, & O'Connell 1993; Saffer et al. 
1994). Webbink (1990, see also Warner 1995) finds a mean 
white dwarf mass of 0.74 ± 0.04 Mq in 84 cataclysmic vari- 
able binaries, and the mean mass of white dwarfs in the 
field is near 0.56 Mq (Bergeron, Saffer, & Liebert 1992). The 
mass of the white dwarf in KPD 0422±5421 (0.62±0.18 Mq) 
is consistent with either of these. These considerations sug- 
gest that the systematic errors are reasonably small and that 
ou derived parameters are trustworthy. 

A close examination of the model light curve near the 
photometric phase 0.5 (when the sdB star is behind the com- 



panion) shows that the white dwarf companion passes in 
front of the sdB star, causing a dip of ss 0.005 magnitudes 
which lasts ps 0.05 in phase (6.5 minutes). There is also a 
very slight distortion at phase 1.0 where the nearly invisible 
white dwarf is totally eclipsed by the sdB star. Interestingly 
enough, the white-light curve from the 0.9m shows a pos- 
sible indication of a small dip near the spectroscopic phase 
0.25 (i.e. the photometric phase 0.5). Caution suggests that 
the observation of a transit be considered only tentative. 
Since the depth of the depression caused by the transit is 
only « 0.005 magnitudes, its exact shape is quite sensitive 
to slight drifts in the transparency, etc. Extended observa- 
tions with a CCD on a large telescope (where one can use 
reasonably short exposures) should be done to confirm the 
existence of the transit. The use of a CCD enables one to 
do precise differential photometry over many orbital cycles. 



5 SUBDWARF ATMOSPHERIC 
PARAMETERS 

The effective temperature T e fj and surface gravity \ogg of 
the subdwarf were determined using a grid of synthetic spec- 
tra. LTE hydrogen and helium line-blanketed model atmo- 
spheres were constructed using version 178 of the atmo- 
sphere code TLUSTY (Hubeny 1988), for solar abundances 
and for reduced helium abundance, [He/H] = —0.5, —1.0, 
and —1.5 (where square brackets denote the logarithmic 
number ratio relative to solar). Detailed spectra were syn- 
thesized on the interval 3780 - 4910 A using version 41 
of SYNSPEC (Hubeny, Lanz & Jeffery 1994); H and He 
I lines were computed using line broadening tables, while 
lines from other elements were computed using Voigt pro- 
files and data from the Kurucz atomic line list. A micro- 
turbulent velocity parameter of 4 km sec -1 was assumed. 
An approximate NLTE treatment of line opacity based on 
second-order escape probability was used in the spectrum 
synthesis (see Hubeny, Harmanec & Stefl 1986). The spectra 
were then convolved and sampled to match the dispersion 
(1.77 A pixel -1 ) and instrumental resolution (3.5 A). On the 
assumption that the subdwarf 's rotation is tidally synchro- 
nized to the binary orbital period, we included a rotational 
broadening of V TO t sini = 150 km sec -1 , which is less than 
the FWHM of the instrumental profile. At each value of 
[He/H], models and spectra were computed at spacings of 
1000 K for T cff = 23000 - 28000 K, and spacings of 0.2 for 
loggr = 5.0 - 5.8. For [He/H] < 0, the grid resolution was 
increased in the neighborhood of the best fit, to 500 K for 
T cft = 24000 - 27000 K and 0.1 for logg = 5.2 - 5.6. 

The models and the restframe spectrum were normal- 
ized to their continua, and the x 2 an d rms residuals were 
computed using most of the spectrum between 3900 A 
and 4900 A. The H/3 profile and nearby continuum were 
not well fit by the models, in part due to the relatively 
poor calibration near the ends. The CII line at 4266.7 A 
is much weaker in the data than in the model spectra. Af- 
ter some experimentation, final fitting regions were chosen 
to be 3949.0 < A < 4210.0, 4282.0 < A < 4394.0, and 
4460.0 < A < 4486.0. This excludes the CII line, some bad 
pixels near 4222 A, and the region near 4430 A which has 
a weak interstellar absorption band. The best-fitting model 
(Figure §) has T cff = 25, 000 K, logg = 5.4, and [He/H] = 
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"Wavelength (A) 

Figure 7. The restframe spectrum (rilled circles) is displayed 
with the best-fitting model (solid line) and the residuals in the 
sense of data minus model (filled triangles), offset by 0.5 units 
upward. The C1I line at 4266.7 A was not included in the fitting 
region. 



— 1.0. The rms residual of the fit is 0.018 (unit = normal- 
ized continuum). For comparison, the value of log<? derived 
from the photometrically determined mass and radius of the 
sdB star is \ogg = 5.50 + 0.21 (Table 5). Keeping T eff , logg, 
and [He/H] at their best values, we find that an approxi- 
mate match of the models to the strength of the CII line 
at 4266.7 A requires a reduction in the carbon abundance 
by a factor of ten or more with respect to solar. Our resolu- 
tion and signal-to-noise ratio are insufficient to make further 
deductions about abundances of metals. 

The estimated sizes of the la statistical errors are 
<tt ~ 1500 K and a"i og£ , ~ 0.1. There is a positive correla- 
tion between the estimated parameters T e a and logg, in the 
sense that forcing a displacement of 0.1 dex in logg requires 
a displacement of ~ 700 K to minimize \ 2 ■ Forcing [He/H] to 
a value away from the best fitting —1.0 pushes T e g and log g 
along this correlation. Our values of T e s = 25, 000 ± 1500 K 
and loggr = 5.4 ± 0.1 are consistent with those values com- 
municated by Saffer (T cff = 26,050, logp = 5.51, [He/H] 
unknown). If we had used SYNSPEC to generate synthetic 
spectra using line-blanketed models from Kurucz (1994 CD- 
ROM distribution) , we estimate that our inferred T e g would 
be ~ 1500 K lower with no change in loggr, based on com- 
parisons made at [He/H]= 0. We stress that we have not 
calibrated our parameter fitting procedure directly against 
stars analysed by other workers. 

The spectroscopic value of logp can be used to make 
an independent estimate of the radius of the sdB star. As a 
function of the total mass the radius of the sdB star is 

1/2 

50 
9 



Rsd 
~llr. 



Mot 
Mr. 



(5) 



M sdB + M comp . On 



where q = M comp /M sdB and A/ tota i 
the other hand, the radius of the sdB star is given by the 
light curve solution as a fraction T of the effective Roche 
lobe radius i?Ri(g). Thus using Kepler's third law, 

1/3 




Total Mass (M ) 

Figure 8. The relations given in Equations ^ and fj are plotted 
as the solid lines. The dotted lines show Equation b using the la 
errors on log g. The dashed lines show Equation using the la 
errors on q The region where the two la regions overlap is outlined 
with the thick line. Finally, the la range of M tota i allowed by the 
mass function is shown with the vertical dashed-dotted lines. 

where Rmiq) is given in Equation ^. In Figure H we plot 
equation q using log g = 5.4 and Equation ^| using q = 0.87. 
The two curves cross at Mtotai ~ 0.7 Mq. We also show the 
lcr limits on each curve (o"i ogg = 0.1 and a q = 0.15). The re- 
gion where the two error bands overlap is outlined with the 
thick lines. An extreme range of 0.25 < Mtotai < 1.75 Ms 
is allowed, compared with the la range of 0.91 < Mtotai < 
1.77 Mq allowed by the mass function and the photomet- 
rically determined values of q and i. We conclude that the 
radius derived from our spectroscopic value of logg is ba- 
sically compatible with the radius derived from the photo- 
metric determination of q and i, while favoring lower values 

Of M to tal- 

A measurement of the projected rotational velocity of 
the sdB star V ro t sin i would provide another measurement 
of its radius: 



R 



Vrot sin i 



RsdB — J'Rni(q) 



P GM U 



47T 2 



(6) 



, 2-7T sin i ) 

provided the sdB star is synchronously rotating. Our present 
spectral resolution is not sufficient to measure Kot sin i for 
KPD 0422+5421 directly. Such a measurement would pro- 
vide another useful constraint on the total system mass 
through the relations given in Equations Fj and pi 



6 THE REDDENING FROM ABSOLUTE 
PHOTOMETRY AND THE DISTANCE 

Absolute photometry of KPD 0422+5421 in the Johnson 
UBV and Stromgren systems was obtained by Downes 
(1986) and Wesemael et al. (1992) respectively. The results 
are: V = 14.66, B -V = +0.20, U - B = -0.65 (no errors 
are quoted); and y = 14.682 ± 0.018, b~y = 0.141 ± 0.011, 
U - b = 0.477 ± 0.001, mi = 0.063 ± 0.019. Since the pho- 
tometric variations are quite small and in phase at different 
wavelengths, it may be assumed that the accuracy of the ab- 
solute photometry is largely unaffected by the cyclical vari- 
ations. 
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Comparison of the U—B and B — V indices with those of 
a large collection of hot subdwarfs given in Fig. 2 of Kilkenny 
et al. (1997b) makes it clear that KPD0422+5421 is anoma- 
lously red. In fact, its colours strongly resemble those of hot, 
subluminous stars with composite spectra, as shown in Fig. 
3 of Kilkenny et al. (1997b). However, our spectra do not ex- 
hibit any signs of anything but the sdB star. Furthermore, 
our photometric solution requires the companion to be de- 
generate. The obvious conclusion is that KPD0422+5421 is 
substantially reddened. 

Further support for this deduction comes from the 
Stromgren colour indices quoted above. Wesemael et al. 
(1992) plot reddening lines in a (u — b, b — y) diagram, 
for a range of assumed sdB temperatures (their Fig. 9); 
KPD0422+5421 inhabits an area in the diagram with E(b — 
y) > 0.2. A crude estimate gives E(b — y) « 0.25, which is 
not surprising given that the star is in the galactic plane. 
By comparing the observed Stromgren photometric indices 
to stellar atmosphere model predictions, Villeneuve et al. 
(1995) derive E(b - y) « 0.26. The value of E{b -y) = 0.25 
implies E(B - V) = 0.36 (Crawford, Glasky, & Perry 1970). 

To compute the distance to the source, we adopt a tem- 
perature of T cff = 25, 000 ± 1500 K, a radius of R sdB = 
0.249+0.008 7?©, a colour excess of E(B-V) = 0.36±0.05, a 
bolometric correction of BC = 2.5±0.1 (Gray 1992), and an 
apparent V magnitude of V = 14.66 (Downes 1986). We find 
Lboi = 21.8+5.4 L Q , M v = 3.86±0.29, and d = 850±130 pc. 
For comparison, Villeneuve, Wesemael, & Fontaine (1995) 
find 581 < d < 1120 pc based on temperatures derived from 
Stromgren photometry. 



7 DISCUSSION OF SIMILAR SYSTEMS AND 
EVOLUTIONARY CONSIDERATIONS 

Three issues relating to the binary evolution arise from our 
study of KPD 0422+5421. 

First, we note that the sdB star, if interpreted as a mem- 
ber of the extreme horizontal branch (EHB), will soon evolve 
to become the second white dwarf in KPD 0422+5421. Thus 
KPD 0422+5421 is a predecessor to the class of double de- 
generate binaries. Presumably it has already passed through 
a common envelope phase. A possible predecessor in turn 
to KPD 0422+5421 is a system like HD 185510 (V1379 
Aql), which is a 21-day binary containing an evolved com- 
pact (sdB) star and a K0 giant (Fekel et al. 1993; Jeffery 
& Simon 1997). The mass ratio of HD 185510 is large, 
Mk/M s( jb ~ 7.5 (Fekel et al. 1993), and one may expect 
that a common envelope stage lies ahead. At that time, HD 
185510 may resemble KPD 0422+5421, with its present sdB 
star in the role of KPD's white dwarf. 

Second, what is the present total binary mass? The 
combined spectroscopic and photometric orbital solution 
suggests a value tantalizingly close to the Chandrasekhar 
limit, Mch ~ 1.4 M©, with all that implies for the possi- 
bility that KPD 0422+5421 may become a Type la super- 
nova, when the components merge. On the other hand, if 
the sdB star is truly a core-He burning EHB star, its mass 
is expected to be close to 0.5 Mq, and with q derived from 
the light curve solution constraining the mass of the unseen 
white dwarf, the total mass may be considerably less than 
Mch- Improved light curves and radial velocity data will 



help to address this question. In the meantime, it is well 
to bear in mind the discussion in Jeffery & Simon (1997) 
concerning whether the canonical 0.5 Mq for EHB stars is 
really empirically well-determined. 

Third, short-period binary systems of the same type as 
KPD 0422+5421, namely sdB + white dwarf, are rare, and 
KPD 0422+5421 arguably possesses the shortest period of 
them all. [One double degenerate system has a shorter pe- 
riod: WD 0957-666, P = 0.060993 days, (Moran, Marsh, 
& Bragaglia 1997). The total mass is estimated to be only 
0.69 M Q , however.] The catalog of Ritter & Kolb (1998) 
lists three possible sdB+ "white dwarf" pairs (i.e. sdB stars 
known to be in close binaries with undetected companions) , 
and the sample of Saffer, Livio, & Yungelson (1998) contains 
seven sdB binaries with undetected companions (including 
two from Ritter & Kolb's catalog), of which five have known 
orbital periods. The six objects with known orbital periods 
(in order of increasing period) are PC 1432+159 (P = 5.39 
hours), PC 2345+318 (P = 5.78 hours), Feige 36 (P = 8.5 
hours), PC 0101+039 (P = 13.7 hours), HZ 22 (=UX CVn, 
P = 13.77 hours), and Ton 245 (P = 2.5 days). A related 
sdO + white dwarf(?) system is HD 49798 with P = 1.55 
days (Thackeray 1970). The best known member of the sdB 
+ white dwarf class is HZ 22 (Humason & Zwicky 1947; 
Young, Nelson, & Mielbrecht 1972; Greenstein 1973) . Green- 
stein (1973) and Young & Wentworth (1982) argue that the 
unseen companion of HZ 22 is a white dwarf. However, it 
is also possible that HZ 22 contains a low mass main se- 
quence star companion, so the nature of HZ 22 is still an 
open question. Another recent sdB+ "white dwarf" candi- 
date is V46 in the globular cluster M4 (Kaluzny, Thompson 
& Krzeminsky 1997), which is an sdB star that shows an ap- 
parently sinusoidal light curve with a period of 1.045 hours. 
This period is much longer than those seen in the pulsating 
sdB stars (EC14026 stars, Kilkenny et al. 1997a). However, 
the nature of the companion star is not known and it is not 
clear if the photometric period is the orbital period. Since it 
has the shortest confirmed orbital period, KPD 0422+5421 
will evolve more quickly via gravitational wave radiation of 
angular momentum than the other sdB + white dwarf sys- 
tems. Ritter (1986) gives the time required for a detached 
binary to reach the semi-detached state: 



t sd = (4.73 x 10yr)- 



M, 



1/3 



,8/3 
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(8) 



where the masses are in solar units, the periods are in days, 
and where 



Psd _ f RsdB \ 

~~P ~ \R m {pS2>) ) • 



(9) 



which assumes that the radius of the sdB star stays fixed. 
Using P sd /P = 0.57 + 0.13, we find that logt sd = 8.17 + 0.15 
(t s d ~ 1.5 x 10 s years), which is comparable to the core He 
burning lifetime of an sdB star (~ 1.5 x 10 8 years, Dorman, 
Rood, & O'Connell 1993). Thus, possibly before the sdB star 
has evolved to a white dwarf, there will be a further episode 
of mass exchange. What KPD 0422+5421 will then look like, 
and whether it can be identified with any presently known 
class of objects, takes the subject well beyond the expertise 
or inclination of the present authors to speculate. 
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8 SUMMARY 

The sdB star KPD 0422+5421 was found to be a short pe- 
riod detached binary with a period of P = 0.0901795 ± (3 x 
10~ r ) days. We argue that the companion star is a white 
dwarf since the companion star is not seen in the spectra, 
and the synthetic light curve model of the U and B light 
curves requires its radius to be on the order of » 0.01 Rq. 
We derive component masses of M s( jb = 0.72 ± 0.26 Mq and 
M comp = 0.62 ± 0.18 M Q . KPD 0422+5421 is one of a small 
number of known sdB + white dwarf binary systems and 
represents a poorly observed and short-lived stage of binary 



star eve Jution 
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